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We use self-consistent field theory to determine structural and energetic properties of intermediates 
and transition states involved in bilayer membrane fusion. In particular, we extend our original 
calculations from those of the standard hemifusion mechanism, which was studied in detail in the 
first paper of this series Q, to consider a possible alternative to it. This mechanism involves non-axial 
stalk expansion, in contrast to the axially symmetric evolution postulated in the classical mechanism. 
Elongation of the initial stalk facilitates the nucleation of holes and leads to destabilization of the 
fusing membranes via the formation of a stalk-hole complex. We study properties of this complex 
in detail, and show how transient leakage during fusion, previously predicted and recently observed 
in experiment, should vary with system architecture and tension. We also show that the barrier to 
fusion in the alternative mechanism is lower than that of the standard mechanism by a few fcsT 
over most of the relevant region of system parameters, so that this alternative mechanism is a viable 
alternative to the standard pathway. 



INTRODUCTION 

The fusion of biological membranes is of great impor- 
tance as it plays a central role inter alia in intracellular 
trafficking, exocytosis, and viral infection 0, E3, 0, H, El ■ 
Given this importance, it might be thought that its mech- 
anism would be well understood, but in fact it is not. 
Perhaps the reason for this is that there is an apparent 
dilemma at the heart of the fusion process. The vesicles, 
or bilayers, which are to be fused must be sufficiently sta- 
ble with respect to irreversible rupture in order to carry 
out their functions on a reasonably long time scale. It 
follows that it must be quite energetically expensive to 
create a long-lived, super-critical, hole in such a mem- 
brane. In other words, the free energy barrier to do so 
must be very large compared to the thermal fluctuation 
energy fcsT. As a consequence, almost all holes created 
by thermal fluctuations do not have sufficient energy to 
traverse this barrier, hence they simply shrink and re- 
seal. However, it is inevitable that for fusion to occur, 
a long-lived hole must be created at some stage of the 
fusion pathway. The dilemma is that a bilayer can both 
be stable with respect to rupture and yet readily undergo 
fusion. 

Some of the solution of this puzzle is in place. It is 
believed that fusion proteins locally expend energy to 
dehydrate both bilayers in order to bring them in close 
proximity. This increases the free energy per unit area 
of the system, i.e puts the system under local stress. As 
a consequence, it is free energetically favorable for the 
system to undergo a transformation that results in a de- 



crease of bilayer area. In principle, this can be accom- 
plished both by fusion and/or rupture, but the proteins 
apparently catalyze the fusion process exclusively. 

The standard hemifusion mechanism, proposed twenty 
years ago by Kozlov and Markin assumes that ther- 
mal fluctuations permit the tails of lipids of the cis leaves, 
those of the apposing membranes which are closest to one 
another, to flip over and form an axially symmetric de- 
fect in the dehydrated region, denoted a stalk [9(. Due 
to the tension, the newly joined cis layers recede so that 
the stalk expands radially preserving the axial symmetry, 
and transforms into a hemifusion diaphragm - a single bi- 
layer consisting of the two remaining trans leaves. Only 
this single bilayer needs to be punctured by a hole in 
order that a fusion pore be formed and the fusion pro- 
cess completed. This radial stalk expansion hypothesis, 
being in qualitative agreement with many experimental 
observations, was essentially the only model of the fusion 
process until recently. 

In contrast to the hemifusion hypothesis, Monte Carlo 
simulations of bilayer fusion 01 showed that fusion can 
evolve through an alternative mechanism [ill Il2| . in 
which the stalk does not expand radially, but rather elon- 
gates in a worm-like fashion. To distinguish the original 
axially symmetric stalk from the elongated structure, we 
will call the former the "classical stalk" for the remainder 
of the paper. Moreover, it was observed that the elon- 
gated stalk destabilizes the fusing membranes by greatly 
enhancing the rate of hole formation in its vicinity. Once 
such a hole is formed in one bilayer close to the elongated 
stalk, the stalk encircles it completely forming a hemi- 
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fusion diaphragm consisting of the other, as yet intact, 
bilayer. Subsequent hole formation in this diaphragm 
completes the fusion process. In a slightly different vari- 
ant of this scenario, holes form in both bilayers near the 
stalk before the stalk has completely surrounded the first 
hole. Fusion is completed when the stalk surrounds both 
holes. This mechanism was also seen in recent Molecular 
Dynamics simulations 0, • ^ was argued that 
the stalk lowers the free energy barrier to hole formation 
by decreasing the effective line tension in that part of the 
hole in contact with the stalk. 

We shall denote the elongated stalk partially surround- 
ing a hole as a stalk-hole complex. As we note below, this 
stalk-hole complex can decay, i.e. evolve without further 
free energy cost, into a final fusion pore, so that this com- 
plex represents a potential transition state in the fusion 
process. 

A direct consequence of this alternative mechanism is 
that there can be transient leakage during fusion. Even 
though leakage is sometimes observed during fusion ex- 
periments, it is usually attributed to the presence of fu- 
sion proteins which are known, for example, to initiate 
erythrocyte hemolysis fl5| . However, the new mecha- 
nism predicts that transient leakage stems from the fu- 
sion pathway itself and should be observable even during 
fusion of model membranes in the absence of fusion pro- 
teins. Leakage during fusion in such systems has indeed 
been observed experimentally 0,0, [3- I n addition, it 
is predicted that this transient leakage should be corre- 
lated in space and time with fusion. Just such correlated 
leakage and fusion were recently observed experimentally 
by Frolov et al. 0|. Fusion without detectable leakage 
is also observed, however El H3- We shall argue 
below that the seeming irregularity of leakage accompa- 
nying fusion can be explained by the new mechanism. In 
particular, the extent of this transient leakage depends 
both on the architecure of the amphiphiles as well as 
the tension (stress) imposed on the membranes. By de- 
creasing the spontaneous curvature of the amphiphiles 
and/or reducing the membrane tension, the leakage can 
be substantially reduced and even completely clliminated 
in some cases. 

While simulations can reveal very clearly the process 
by which fusion occurs, it is very difficult to extract free 
energy barriers involved in the process. Moreover, it 
would be prohibitively expensive to study comprehen- 
sively the fusion process in a range of parameters such as 
amphiphile architecture and/or bilayer tension. 

In a previous paper Q, we employed self-consistent 
field theory (SCFT) to evaluate these barriers assum- 
ing that fusion took place via the standard, radially- 
expanding stalk and hemifusion mechanism. The system 
considered consisted of bilayers of AB block copolymer, 
with fraction / of the A monomer, in a solvent of A 
homopolymer. All polymers were characterized by the 
same polymerization index. Comparison of various prop- 



erties of this system of block copolymer amphiphiles with 
those of membranes consisting of biological lipids permit- 
ted an estimate that free energies of a structure in the 
copolymer simulations were 2.5 times smaller than those 
of the corresponding structure in the biological system. 
We calculated the barrier to stalk formation in polymeric 
bilayers, and from it estimated that in membranes made 
of biological lipids, this barrier would not exceed 13fcsT. 
The larger barrier in the standard process is that asso- 
ciated with the radial expansion of the hemifusion di- 
aphragm [2^, and we estimated this to be in the range 
of 25-63fcsT, depending upon the lipid architecture and 
membrane tension. Perhaps one of the most interesting 
results of this study was the following: the range of varia- 
tion in amphiphile architecture over which successful fu- 
sion can occur is severely restricted by the fact that the 
fusion process begins with the formation of a metastable, 
classical stalk. If / is too large, corresponding to lipids 
with very small spontaneous curvature, stalks between 
bilayers are never metastable. On the other hand, if / 
is too small, corresponding to lipids with larger nega- 
tive spontaneous curvatures, linear elongated stalks be- 
came favorable which destabilize the bilayers completely 
by causing a transition to an inverted hexagonal phase. 
Thus in order for fusion to occur, the lipid composition 
of membranes must be tightly regulated. This conclusion 
also applies to fusion which proceeds via the new mecha- 
nism as it, too, begins with the formation of the classical 
stalk. 

In this paper we apply SCFT methods to calculate the 
fusion barriers in this new, alternative, mechanism. We 
begin in Section D by calculating the free energy of an iso- 
lated hole in a single bilayer as a function of its radius, 
R, for bilayers under various tensions, 7 and consisting 
of diblocks of different architectural parameters, /. The 
result is that, as expected, it is very expensive to create 
a hole in an isolated bilayer. In Section B we turn to the 
calculation of the free energy of the stalk-hole complex. 
Because this complex is not axially symmetric, our task is 
much more difficult than our previous calculation of the 
barriers in the old hemifusion mechanism in which the 
intermediates were postulated to be axially symmetric. 
We accomplish our goal by constructing the non-axially 
symmetric intermediates from fragments of other excita- 
tions which do possess this symmetry, and therefore are 
more easily obtained. In particular, we show that struc- 
tures related to an Inverted Micellar Intermediate (IMI) 
play important role in this process. We compare our re- 
sults for the free energy barrier in the two mechanisms 
and show that the barrier in the new one is indeed lower 
than that in the old, although the difference in most of 
the region of parameters in which fusion can occur suc- 
cessfully is not more than a few fe^T. Finally in Section B 
we discuss further the reason why the new mechanism 
is a favorable one. We trace it not only to the reduc- 
tion of the line tension of a hole when nucleated next to 
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a stalk, but also to the relatively low cost for the stalk 
to extend linearly. Consequently when a hole appears 
in the bilayer, a large fraction of its circumference can 
have its line tension reduced by the nearby presence of 
a stalk. We conclude with some comments on the de- 
pendence of the rate of hole formation in a bilayer on 
the line tension of the hole. We show that even modest 
changes in the effective line tension of a hole due to the 
presence of the elongated stalk in the stalk-hole complex 
can strongly affect the rate of hole formation, and hence 
the rate of fusion. Such small changes in line tension, 
therefore, destabilize what were very stable bilayers and 
enable them to undergo fusion. 

FREE ENERGY OF A HOLE IN AN ISOLATED 
BILAYER 

In this section we discuss the free energy of a circu- 
lar hole in an isolated bilayer to show that the energy 
associated with formation of such a defect is high, as is 
expected if isolated bilayers are stable. The SCFT cal- 
culation follows the lines described in our previous paper 
1] . It is straightforward within the SCFT to obtain the 
free energy, fi m (T, A/x, V, A), of a bilayer of area A at a, 
temperature T and a difference, A/x = /x Q — /x s , of the 
bulk chemical potentials of the amphiphile and of the 
solvent. There is only one independent chemical poten- 
tial as the system is assumed to be incompressible. The 
volume of the system is V. Similarly, we denote the free 
energy of the system without the bilayer, i.e., a homoge- 
neous amphiphile solution, fla(T, A/x, V). The difference 
between these two free energies, in the thermodynamic 
limit of infinite volume, defines the excess free energy of 
the bilayer membrane: 

dn m (T,Afi,A)= lim [n m (T,An,V,A)-n (T,A[x,V)}. 

V—>oo 

(1) 

The excess free energy per unit area, in the thermody- 
namic limit of infinite area, defines the lateral membrane 
tension 

7 (T,A/x)= lim [8n m (T,Afx,A)/A]. (2) 

A— >oo 

Changes in this tension 7 can be related to changes in 
the temperature and chemical potential by means of the 
Gibbs-Duhcm equation 

d 7 (T, A/x) = -Ss AT - Sa a d(Afi), (3) 

where 5s is the excess entropy per unit area, and 5a a 
is the excess number of amphiphilic molecules per unit 
area. This relation shows that the chemical potential 
difference, A/x, can be used to adjust the bilayer tension 
7- 

As discussed previously it is also possible to in- 
troduce axially symmetric defects of a specified radius R 
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FIG. 1: Density profiles of bilayers pierced by an isolated hole 
are shown for three different hole radii: R/R g = 1, 2, and 5. 
Only the majority component is shown at each point. Solvent 
segments are white. Hydrophilic and hydrophobic segments 
of the amphiphile are shaded dark and light respectively. 



into the bilayer and to obtain the excess free energy of 
such structures. The choice of model parameters was dic- 
tated by our original Monte Carlo simulations 01 an d 
the details can be found in the first paper of this series 


Fig. shows the density distribution of hydrophobic 
(B) and hydrophilic (A) segments in a bilayer with holes 
of different radii, which are defined as the radial distance 
in the plane of symmetry to the A/ B interface, the point 
at which the volume fractions of A and B monomers are 
equal. We find that qualitative features of this profile 
are not very sensitive to the architectural parameter / or 
tension 7: the rim of the hole has a shape of a bulb which 
is typical whenever a flat bilayer has an edge |2~1. 125L Etjf . 

The free energy of such a hole in a bilayer is shown 
in Fig. EJa) as a function of its radius for a bilayer at 
zero tension and composed of amphiphiles of different 
architectural parameters, /. One sees that under zero 
tension the free energy increases essentially linearly with 
R and the excess free energy of the hole can be writ- 
ten as 2ttXh(T,Aii,R)R, where Xh is an effective line 
tension. As one would expect, this line tension quickly 
asymptotes to a constant value A#(T, A/x) for sufficiently 
large R. For the bilayer under zero tension composed of 
amphiphiles with / = 0.35, we find XnRg/kT = 2.63. 
To compare with analogous values for membranes, we 
convert this to the dimensionless ratio Xu/lod, where 
70 is the free energy per unit area of an interface be- 
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FIG. 2: (a) Free energy of a hole in an isolated bilayer as 
a function of R/R g at zero tension for various amphiphile 
architectures, /. From top to bottom the values of / are 0.29, 
0.31, 0.33, and 0.35. (b) Same as above, but at fixed / = 0.35 
and various tensions 7/70- From top to bottom, 7/70 varies 
from 0.0 to 0.6 in increments of 0.1. 

tween coexisting phases of bulk homopolymer A and bulk 
homopolymer B, and d is the thickness of the bilayer. 
From our previous work jioj . we obtain 7oc£ 2 /kT = 65.3, 
and d/Rg = 4.47 so that X H /j d = 0.18. The analo- 
gous quantity can be calculated for membranes taking 
A™ = 2.6 x lCT 6 dynes |H d m = 35.9 x 1CT 8 cm 
[29]. and an oil- water tension of j m = 50 dynes/cm, from 
which X m /"fmd m = 0.14. Thus the line tensions we ob- 
tain are reasonable. 

In Fig. |2tb) we show the effect of membrane tension 
on the dependence of the hole free energy on radius R. 
One sees, as expected, that the free energy of the hole 
eventually decreases due to the elimination of stressed 
membrane area. For sufficiently large R, one expects the 
free energy of the hole to be of the form 

F H (R) = 2tt\ h R- 7r 7 i? 2 , (4) 

with 7 the imposed membrane tension. We verified that 
this form is certainly adequate at large R. At smaller 
radii, which will be of interest to us, the coefficients Xh 
and 7 are, themselves, functions of R. One sees from this 
figure that under relevant tensions the maximum value 
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FIG. 3: Parametrization of the elongated stalk. Gray 
schematically shows location of the hydrophobic segments in 
the plane of symmetry between fusing bilayers. The arc ra- 
dius R corresponds to the radial distance to the outer hy- 
drophilic/hydrophobic interface in the plane of symmetry. 
Values of the fractional arc angle a, defined in the range [0, 1], 
are given at the top of each stalk configuration. Note that 
a = corresponds to the original stalk, whereas a = 1 cor- 
responds to a family of structures reminiscent of the IMI (see 
also Fig. |3). 

of Fh(R), which is required to form a critical size hole 
leading to irreversible membrane rupture, is no less than 
about 16k bT. In a bilayer composed of biological lipids, 
we can estimete that this would correspond to a barrier 
of approximately AOksT. Thus, as stated, isolated bi- 
layers are very robust against rupture caused by thermal 
excitation, and it is precisely this stability that makes 
fusion difficult to understand. 

We now turn to the calculation of the stalk-hole com- 
plex, which is a possible fusion intermediate. We will 
show that the barrier to fusion is much less than the bar- 
rier to create a hole in each of the two fusing bilayers in 
the absence of an elongated stalk. 



FREE ENERGY OF THE STALK-HOLE 
COMPLEX 

Immediately before formation of stalk-hole complex 

Right after the formation of the initial (classical circu- 
lar) stalk and just before the formation of the stalk-hole 
complex, the stalk elongates in a worm-like fashion. For 
the sake of simplicity, we assume that in the z — (sym- 
metry) plane this elongated structure has a shape of a 
circular arc with a fractional angle, < a < 1, and 
radius R, as shown schematically in Fig. With this 
choice of the parameters, a = corresponds to the clas- 
sical stalk structure, whereas for a = 1 there is a family 
of structures that are reminiscent of the Inverted Micel- 
lar Intermediate (IMI) , studied previously by Siegel [3(| • 
Although the radii of the structures we consider are not 
necessarily small, we shall continue to refer to them as 
IMIs. A density profile of one such structure is shown in 
Fig.0| Its radius R is defined as the radial distance to the 
furthcrest point on the z — plane at which the densi- 
ties of hydrophobic and hydrophilic segments are equal, 
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FIG. 4: Density profile of an IMI. The amphiphiles are char- 
acterized by / = 0.3. The radius of the IMI, in units of the 
ideal radius of gyration, R/R g is 3.4. Gray scale as in Fig.Q 



and is shown in the figure. We denote its free energy 
Fjmi(R)- Note that the equilibrium IMIs considered by 
Siegel correspond to structures with an optimal radius 
R* , which minimizes Fimi(R). 

In general, the elongated stalk will not form a complete 
IMI, that is a will be less than unity, so we approximate 
the free energy of the extended stalk in this configuration 
as 



F 1 (R,a) = aF IM i{R)+F s . 



(5) 



The presence of the second term is due to the free energy 
of the end caps of the extended stalk (see Fig. I3J). As 
these two ends together form an axially symmetric stalk, 
the free energy of these ends is just the free energy of the 
classical stalk, F s , which we have calculated previously 
0. Note that for the case a = 1, the above estimate is 
certainly an upper bound as the second term should be 
absent in that case. 

The free energy of the IMI can be calculated readily 
because it possesses axial symmety. In Fig.^a) we show 
its free energy as a function of radius for a bilayer un- 
der zero tension for various architectural parameters /. 
Again, as is the case with the other axially symmetric 
structures we studied, the free energy is asymptotically 
linear at large R, with the slope 2-kXimi defining the ef- 
fective line tension Xjmi ( c -/- Eq. QJ. In Fig. Efb) we 
also show the free energy of the IMI as a function of R 
for a bilayer with fixed / = 0.31 and different tensions. 
From[3]it is apparent that the free energy of the structure 
for the sizes that are pertinent to the fusion intermediate 
cannot be described by a simple estimate based on the 
line tension of the IMI. The increase of the free energy 
with decreasing radius at small radius results from the 
repulsion of the interfaces across the IMI structure. It 
is similar to the free energy barrier associated with clos- 
ing the fusion pore 0. Note that the free energy does 
not decrease with R for large R because the IMI does 
not eliminate bilayer area. Therefore for large enough a 
and/or radius R, the free energy of this structure will 




FIG. 5: (a) Free energy of an IMI as a function of R/R g at 
zero tension for various amphiphile architectures, /. (b) Free 
energy of an IMI with / = 0.31 and for various tensions 7/70- 
From top to bottom, 7/70 varies from 0.0 to 0.4 in increments 
of 0.1. Minima on these curves correspond to metastable IMI 
structures. 



exceed that of the stalk-hole complex in which a hole 
forms next to the elongated stalk. We turn now to the 
calculation of the free energy of this complex. 



Immediately after formation of stalk-hole complex 

We model the stalk-hole complex as an elongated stalk 
in contact with a circular hole in one of the bilayers. 
We assume that the radial axes of the elongated stalk 
and of the hole coincide and that the radius of the hole 
is R — 5, where R is the radius of the elongated stalk, 
and that S is chosen such that the hole is aligned in the 
radial direction with the elongated stalk over a fraction 
of its circumference, again denoted by a. See Fig. HjJ To 
calculate the free energy of this configuration, we note 
that at a = 1, the configuration is simply a hcmifusion 
intermediate (HI) of radius R, and the elongated stalk 
would now connect two bilayers to one bilayer. We have 
calculated the free energy of the hemifusion intermediate 
previously 0. The radius R of this structure is defined 
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FIG. 6: Parametrization of the stalk-hole complex. Gray 
schematically shows location of the hydrophobic segments 
in the plane of symmetry between fusing bilayers. The 
arc radius R corresponds to the radial distance to the hy- 
drophilic/hydrophobic interface of the hemifusion intermedi- 
ate in the plane of symmetry. Projection of the edge of a hole 
in one of the membranes is shown with dashed line. The ra- 
dius of this hole is R — 8. The other membrane does not have 
a hole. The hydrophobic thickness of the bilayer is 5. Values 
of the fractional arc angle a, defined in the range [0,1], are 
given at the top of each stalk configuration. 



by the position of the hydrophilic/hydrophobic interface 
in the z = symmetry plane. With these definitions of 
the radii of the hemifusion intermediate and of the hole, 
a choice of 6("f, /) equal to the hydrophobic thickness of a 
bilayer, ensures that the hole is adjacent to the elongated 
stalk. In general, when the hole forms, the elongated 
stalk does not completely surround it, so that a fraction 
1 — a of the stalk-hole complex looks like a bare hole 
edge in an isolated bilayer. Thus we approximate the 
free energy of this stalk-hole complex to be 

F 2 (R,a)=aF HI (R) + (l-a)F H (R-S)+F d . (6) 

The free energy F d comes from the end caps of the elon- 
gated stalk connecting to the hole edge. The two ends 
together do not make an axially symmetric stalk, but like 
the stalk, this defect is also saddle-shaped so one expects 
its free energy to be small and not very different from 
that of the stalk. 



The transition state 

It is clear from Eqs. [5] and that the free energies of 
these structures depend both on the radius, R, of the 
intermediate and on the fraction, a. Thus we must con- 
sider a two-dimensional reaction coordinate space, (R, a). 
The fusion process starts off by the formation of the clas- 
sical stalk, which corresponds to the a — line on the 
Fi(R, a) free energy surface. Elongation of the stalk cor- 
responds to non-zero values of a. We assume that the 
stalk-hole complex forms when the free energy surfaces 
Fi(R,a) and F 2 (R,a) intersect. This intersection hap- 
pens along a line in the (R, a) plane, which defines the 




FIG. 7: Four free energy landscapes (in units of ksT) of 
the fusion process, plotted as a function of the radius, R (in 
units of R g ) and circumference fraction a. The architecture of 
the amphiphiles and the value of the tension 7/70 are given. 
The dotted line shows a ridge of possible transition states, 
separating two valleys. The region close to the a = line cor- 
responds to a barely elongated stalk intermediate (see Eq.^jJ. 
The other valley, close to a = 1 states, corresponds to a hole 
almost completely surrounded by an elongated stalk. The sad- 
dle point on the ridge, denoted by white dot, corresponds to 
the optimal (lowest free energy) transition state. The energy 
of the defect, Fd has been set to zero here. 



"ridge" of possible transition states (R, aTs(R)) with 



a T s(R) 



= 1- 



F H (R-S)+F d -F s 
F H (Rs) + F IM i(R) - F H i(R) ' 

Fimi{R) — Fhj(R) + Fs — Fd 



F H (R -5) + F IMI (R) - F HI {R) 



(7) 



(8) 



The free energy of the optimal transition state can ob- 
tained by finding the free energy minimum along the 
ridge of the transition states, which we shall do momen- 
tarily. First we note from Eq. [8] that the fraction of the 
hole surrounded by the elongated stalk increases as the 
free energy of an isolated hole increases. This shows that 
the reduction of the high cost of the bare hole edge is a 
driving force of this mechanism. 

We return to the free energy landscape of the fusion 
process defined by min(Fi(i?, a),F 2 (R, a)). Examples of 
such landscapes are shown in Fig. To clarify the effect 
of different parameters we present results for a membrane 
consisting of lipids with / = 0.31 and 0.33, and under the 
reduced membrane tension 7/70 = 0.1 and 0.2. To ob- 
tain these results we have set the small defect energy F d 
to zero. This parameter has very little effect on the qual- 
itative features of the landscapes. Quantitative effects 
are also small and will be discussed below. 

The landscapes are saddle-shaped with low free energy 
valleys close to a = and a = 1 lines. The first valley 
corresponds to barely elongated stalks of very small cir- 




FIG. 8: Plot of a* , which corresponds to the optimal transi- 
tion state in the stalk-hole mechanism, as a function of archi- 
tecture of the amphiphiles and the tension of the membrane. 



cumfcrence, configurations which are clearly energetically 
inexpensive. The second valley corresponds to a hole that 
is almost completely surrounded by the elongated stalk. 
Its energy is small because formation of the hole leads 
to a decrease of the membrane area under tension. One 
should note that a = 1 corresponds to the hcmifusion 
intermediate, which is also formed in the standard mech- 
anism, but through a completely different pathway. 

The ridge of the transition states (R,o.ts(R)) is indi- 
cated by a dotted line. There is a saddle point along this 
ridge, denoted by a circle on the plots. We denote the 
value of the radius of this optimal transition state as R* , 
and the value of ceTs{R*) as a*. The free energy of the 
transition state F* = Fx(R*,a*) = F 2 (R*,a*) This as- 
sumes that one can ignore any additional barriers caused 
by the rearrangement of amphiphiles in passing from the 
configuration just before the formation of the stalk-hole 
complex to the configuration just after. 

The value of a at the saddle point, a*, is shown in 
Fig. |SJ Once the stalk-hole complex has formed, the free 
energy of the complex decreases as the stalk continues to 
enclose the hole, that is, as a increases to unity. This is 
clear a priori because as the stalk advances around the 
perimeter of the hole, it reduces the large line tension 
of the bare hole to the smaller line tension of the hole 
surrounded by stalk without any concommitent increase 
in energy due to that advance. 

For small values of the architectural parameter /, there 
is a considerable region for which a* — I. The reason for 
this can be inferred from Fig. [5] which shows the calcu- 
lated asymptotic (large R) values of the line tension, Xes, 
of the elongated stalk. One sees that Xes decreases as a 
function of / so that the free energy of an IMI, Fjmi(R), 
which is dominated by this line tension, also decreases. 
Thus when the hole appears next to the extended stalk in 
a membrane characterized by a small /, more of the hole 
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FIG. 9: Line tensions of an elongated linear stalk, Xes, of a 
bare hole in a membrane, \h , and of a hole which forms next 
to an elongated stalk, Xsh as a function of architecture, /. 
All line tensions are in units of ksT / R g . 

will be surrounded by the stalk, that is, a will increase 
towards unity. This physical explanation is reflected in 
Eq.El 

We expect that this result has consequences for the 
amount of transient leakage during the fusion event. It 
is reasonable to expect that the amount of leakage would 
decrease as (1 — a), because this is the fraction of the hole 
in the stalk-hole complex which is not "sealed" by the 
stalk. In fact, for architectures with sufficiently small /, 
(i.e. sufficiently large, negative spontaneous curvatures), 
Fig. leads us to expect that (1 — a) = so that there 
would be no transient leakage at all. 

The free energy barrier to formation of the stalk-hole 
complex measured relative to the initial metastable stalk, 
(F* — Fs)/Js;bT, is shown in Fig. ITUT a). For comparison, 
we also show the barrier encountered in the standard 
hcmifusion expansion mechanism, which we calculated 
earlier for the same parameters. It is clear that in both 
mechanisms the free energy barrier can be significantly 
lowered either by an increase in the membrane tension 
or decrease in the hydrophylic fraction / (more negative 
spontaneous curvature). The difference between these 
two barrier heights, in units of fc^T, is shown in Fig. ITTI 
It is positive when the barrier in the old mechanism ex- 
ceeds that of the new mechanism. We see that over the 
entire region, the barrier to fusion is lower in the new 
mechanism, and becomes increasingly favorable as / de- 
creases, i.e. as the amphiphile architecture becomes more 
inverted-hexagonal forming. We estimate that the differ- 
ence in barrier heights in this system of block copolymers, 
from about 1 to 7ksT, would translate to a range of 3 to 
18/cbT in a system of biological lipids. 

As noted earlier, we have set Fd to zero in the above. 
Recall that the defect is the free energy of the two caps 
at the ends of the rim of an incomplete hemifusion di- 
aphragm. As these caps are similar, but not identical, to 
the two halves of a stalk, we expect their energies to be 
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FIG. 10: Free energy barriers, in units of fc_gT, in (a) the new 
stalk-hole complex mechanism, and (b) the standard hemifu- 
sion mechanism. 



similar. From our calculations 0, we know that the stalk 
free energy in our system does not exceed 4fcsT. If we set 
Fd to 4fc gT, then the difference in barrier heights in the 
two mechanism changes somewhat, and is shown in part 
(b) of Fig. ^] The new mechanism is still favored over 
most of the tension/architecture space, while the stan- 
dard mechanism is now favored for bilayers composed of 
amphiphiles with larger values of / i.e lamellar-forming 
lipids and under small tensions. 



Formation of the final state 




0.35 



0.35 



FIG. 11: (a) Difference between the free energy barrier in the 
standard mechanism and that in the new mechanism, in units 
of fesT, as a function of architecture and tension. The defect 
free energy is here taken to be zero, (b) Same as in (a) except 
that the defect free energy is taken to be 4fc bT. 



small region at small / and low tension, shown in black 
in Fig. El This region of parameters corresponds to the 
stalk-hole transition state with a = 1, i.e. a completely 
formed IMI (see Fig.[5J). Therefore, we conclude that for- 
mation of the stalk-hole complex involves the largest free 
energy barrier along this pathway. In the special situa- 
tion of very negative spontaneous curvature amphiphiles 
and low tension, the system does not continue on to the 
formation of a fusion pore, but remains in a state in which 
the membranes are joined by an IMI structure. 



The stalk-hole complex is a transition state along the 
fusion pathway, but for complete fusion to occur it has 
to transforn into the fusion pore. Properties of the fusion 
pore have been considered in detail in our first paper Q . 
In the case of the hemifusion mechanism we have found 
that the fusion pore has a lower free energy than the 
hcmifused transition state and thus can be formed with- 
out an appreciable additional barrier. In the present case 
we also found that the fusion pore has a lower free en- 
ergy than the transition state everywhere except a very 



DISCUSSION 

We have utilized self consistent field theory and a 
model of polymeric bilayers to calculate the free energy 
barriers along the fusion pathway first seen by Noguchi 
and Takasu and by ourselves [ill fl2j . There are at least 
two barriers associated with this path; a smaller one asso- 
ciated with the formation of the initial axially symmetric 
classical stalk, and a larger one associated with the for- 
mation of the stalk-hole complex. This path replaces the 
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FIG. 12: Difference in free energy, in units of fcsT, between 
the stalk-hole transition state and fusion pore of the same 
radius. 

expensive step in the old mechanism, which is the radial 
expansion of the stalk into a hemifusion diaphragm, by 
the expensive step of elongating the stalk in a worm-like 
fashion and having a hole form next to it which creates 
the stalk-hole complex. There are several points that we 
wish to make. 

First, by direct comparison of the calculated free en- 
ergy barriers in the new mechanism and in the standard 
one, we have demonstrated that the free energy barriers 
are comparable. Hence this new pathway is a viable alter- 
native to the standard mechanism. We have also demon- 
strated the tendency that the new mechanism tends to be 
the more favorable the more the amphiphile architecture 
approaches that of inverted-hexagonal formers. 

Second, as noted previously, the new mechanism pre- 
dicts the possibility of transient leakage which is corre- 
lated in space and time with fusion. Just such leakage, 
correlated in space and time with fusion, has been ob- 
served 0| . This prediction is in contrast with the old 
mechanism in which any leakage that occurs is not corre- 
lated directly with the fusion process itself. Our calcula- 
tions predict that the amount of this correlated leakage 
decreases, and can vanish altogether, as the architecture 
of the amphiphiles becomes more like that of inverted- 
hexagonal formers. This is a prediction that could be 
tested by carrying out a series of experiments like those 
of Frolov et al. |l9j on vesicles for which one could vary 
the amphiphile architecture or the relative composition 
of amphiphiles of different architecture. Such control of 
amphiphile architecture is readily obtained in polymer- 
somes ,31j which would therefore offer an excellent system 
in which to test this prediction. 

Third, our calculations predict existence of equilibrium 
IMI structures that are metastable with respect to for- 
mation of the fusion pore, except in a region of very low 
/ and 7 where they are actually favored over a pore. 
The possible occurrence of these structures had previ- 



ously been dismissed due to very high estimates of the 
free energy of their formation |30j . 

Finally we observe that in order for this new mecha- 
nism to be favorable, two conditions must be met. The 
first is that it must not cost too much free energy for the 
stalk to elongate in a worm-like fashion, in the manner 
that it does before the hole appears. That this can be 
the case is clear from the fact that at the transition to an 
inverted hexagonal phase, the line tension of linear stalks 
is small. Thus as the architecture is varied such that the 
system approaches this transition, it must be inexpensive 
for the stalk to elongate and wander. That this is correct 
can be seen from the calculated line tension, XeS: 01 the 
elongated linear stalk shown in Fig. El It is essentially 
independent of tension, 7. We see that this line tension 
decreases with decreasing / as expected, which decreases 
the cost of elongating a stalk. The second condition is 
that the free energy of the hole which is created must 
not be too large. As noted earlier, the high cost of an 
isolated hole is due to the line tension of its periphery. 
If this is reduced by causing the hole to form next to 
the elongated stalk, the cost of the hole in the stalk-hole 
complex will also be reduced. To determine whether this 
is so, we have calculated the line tension of an isolated 
hole in a bilayer, Ah, and also the line tension of a hole 
created next to an elongated stalk, Ash- These results, 
again essentially independent of the membrane tension, 
are shown in Fig. |5J as a function of architecture. It is 
seen that in the region of / in which successful fusion is 
possible, 0.29 < / < 0.37 [l|, the line tension of the hole 
is reduced by about a factor of two. Let us now show 
that even such a relatively small change can have a very 
large effect on the rate of fusion. 

Consider the simple estimate of the free energy of a 
hole, Eq. which we reproduce here 

F H = 2n\ H R-ir-fR 2 . (9) 

The height of the barrier to stable hole formation cor- 
responds to the maximum of this function. We ignore 
any i?-dependence of Ah and 7 and immediately ob- 
tain the radius of the hole corresponding to the barrier 
to be R* = Aff/7, and the height of the barrier to be 
F* = ttA h /j. The rate of formation of an isolated hole 
in a bilayer is proportional to the Boltzmann factor 

P H = e^{-[F* -k B T\n(A H /i 2 )]/k B T}, (10) 
= if expt-vrA^MsT), (11) 

where the entropy associated with the formation of a hole 
in an available area Ah is — fcs hi(An /t 2 ) with i a char- 
acteristic length on the order of the bilayer width. If 
Ph <C 1, then the bilayer is stable to hole formation by 
thermal excitation. 

The formation of the stalk-hole complex reduces the 
line tension of that part of the hole near the stalk from 
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\h to Xsh- This can be described by introducing the 
effective average line tension entering Eq. 1111 



A 



H 



Xr 



a\sH + (1 - u)X H - 



(12) 



Then the corresponding rate of stalk-hole complex for- 
mation becomes 



Pi 



SH 



Nsas 
i 2 



exp (-nXl hk B T) , 



(13) 



where N$ is the number of stalks formed in the sys- 
tem and as is the area around each stalk in which 
hole nuclcation can take place. For the small reduction 
Xsh/Xh = 1/2, the above becomes 



P.' 



SH 



Ph 



Nsas 
A H 

N s as 



exp 



"fk B T 



A \ a(l-a/4) 



e 2 p, 



H 



> (14) 
(15) 



This shows explicitly that if the isolated membrane is 
stable to hole formation, (i.e. Ph <C 1), then even a small 
reduction in the line tension ensures that formation of the 
stalk/hole complex causes the rate of hole formation in 
the apposed bilayers, and therefore fusion, to increase 
greatly. 

We illustrate this with two examples. We first consider 
the copolymer membranes which we simulated previously 
In that case the exponent in the Boltzmann 



tne cop < 
factor 



jk B T 



XffRg 

k H T 



7o\ (k B T\ 



(16) 



where 70 is the tension of an interface between bulk hy- 
drophilic and hydrophobic hompolymer phases. The var- 
ious factors in the simulated system are X B Rg/k B T = 2.6 
at / = 0.35 (see FigJ^Land 7o/7 = 4/3, k B T/ l0 R 2 g = 



0.31, A H /£ 2 = 39 |10l Il2j. From this we obtain P H 
6 x 10~ 3 , so that isolated bilayers should have been sta- 
ble to hole formation, as was indeed the case. However 
in the presence of a stalk, the Boltzmann factor will 
be increased according to Eq. If we assume that 

the elongated stalk enclosed one half of the perimeter 
of the hole when it appeared, (i.e. a — 1/2), and that 
nsas /A B ~ 0.3 (consistent with the simultaneous obser- 
vation of multiple stalks in a small simulation cell |l2j). 
we find that Psh/Ph ~ 14 so that the rate of hole for- 
mation should have increased appreciably as observed in 
the simulations. 

This increase is expected to be more dramatic in bio- 
logical membranes. In that case we estimate the expo- 
nent of the Boltzmann factor, —irX H /jk B T, as follows. 
We take the line tension to be that measured in a stearoy- 
loleoylphosphatidylcholine and cholesterol bilayer, X B ~ 
2.6 x 10 -6 erg/cm [2^,|2^. For the surface tension, we 



take an estimate of the energy released by the conforma- 
tional change of four of perhaps six hemagglutinin trimers 
aranged around an area of radius 4nm, each trimer giv- 
ing out about 60k B T [32j . This yields an energy per unit 
area 7 w 20 erg/cm 2 . Thus P H = 1.7 x 10- 11 (A H /£ 2 ), 
which indicates that even subject to this large, local, en- 
ergy per unit area, the membrane is quite stable to hole 
formation for vesicles of any reasonable size. However if 
we assume again that the line tension of the hole is re- 
duced by a factor of two by being nucleated next to the 
elongated stalk, that the stalk extends halfway around 
the circumference of the hole, and the density of stalks is 
such that nsas /A B = 0.3, then the rate of hole forma- 
tion is increased by 



P S H 
Pi 



= 0.3 



H 



1 



1.7 x 10" 



7/16 



1 x 10 4 



(17) 



or an increase of more than four orders of magnitude. 

One should note the implications of this simple argu- 
ment. Because the probability to form a stable hole de- 
pends exponentially on the square of the line tension, an 
isolated bilayer is guaranteed to be stable against hole 
formation for normal line tensions. However it is pre- 
cisely this same dependence which also ensures that the 
bilayer will be destabilized by hole formation due to any 
mechanism which even modestly reduces that line tension. 
From here it is only a short step to successful fusion. 
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